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ABSTRACT 

We present probability distribution functions (PDFs) of the surface densities of 
ionized and neutral gas in the nearby spiral galaxies M31 and M51, as well as of dust 
emission and extinction Ay in M31. The PDFs are close to lognormal and those for 
Hi and Ay in M31 are nearly identical. However, the PDFs for H 2 are wider than 
the Hi PDFs and the M51 PDFs have larger dispersions than those for M31. We use 
a simple model to determine how the PDFs are changed by variations in the line-of- 
sight (LOS) pathlength L through the gas, telescope resolution and the volume filling 
factor of the gas, fy. In each of these cases the dispersion a of the lognormal PDF 
depends on the variable with a negative power law. We also derive PDFs of mean 
LOS volume densities of gas components in M31 and M51. Combining these with the 
volume density PDFs for different components of the ISM in the Milky Way (MW), 
we find that cr decreases with increasing length L with an exponent of —0.76 ± 0.06, 
which is steeper than expected. We show that the difference is due to variations in 
fy. As fy is similar in M31, M51 and the MW, the density structure in the gas in 
these galaxies must be similar. Finally, we demonstrate that an increase in /„ with 
increasing distance to the Galactic plane explains the decrease in cr with latitude of 
the PDFs of emission measure and FUV emission observed for the MW. 

Key words: ISM: structure - turbulence. Galaxies: individual - M31, M51, MW 


1 INTRODUCTION 

Many processes can leave an imprint on the complicated 
structure of the interstellar medium (ISM) observed in the 
Milky Way and nearby galaxies. For example, variations in 
the gas density can be due to: an evolving population of 
gravitationally bound clouds; expanding supernova remnant 
shells compressing the gas; cooling and heating of the gas 
leading to phase transitions; the turbulent and compressible 
nature of the flow. 

During the last two decades compressible turbulence 
has been recognized as an important cause of the density 
structure of the ISM in the Milky Way. Both isothermal 
( Elmegreen Sz Scalo|2004 and references therein) and multi¬ 
phase (de Avillez & Breitschwerdt 2005| Wada & Norman 
|2007[ ) simulations of the ISM produced lognormal volume 
density distributions. Several authors noted that for a large 
enough volume and for a long enough simulation run, the 
physical processes causing the density fluctuations in the 
ISM in a galactic disc can be regarded as random, inde¬ 
pendent, multiplicative events (e.g. Vazquez-Semadenijl994 


[Passot fc Vazquez-Semadeni|[T998[ |Wada fc Norman||2007p . 


* E-mail: eberkhuijsen@mpifr-bonn.mpg.de 
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Therefore, the PDF of log(density) becomes Gaussian and 
the density PDF lognormal, as long as the local Mach num¬ 
ber and the density are not correlated ([Kritsuk et alT]|2007[ 


Federrath et al. 20101. However, Tassis et al. (20101 show 


that lognormal column density distributions in molecular 
clouds are a generic statistical representation of the inho¬ 
mogeneous distribution of gas density, whether the inhomo¬ 
geneity is a result of turbulence, gravity or ambipolar diffu¬ 
sion. Recent reviews which discuss density PDFs of molec¬ 
ular clouds and their connection to turbulence and star for¬ 
mation includelMcKee & Ostriker (20071; Hennebelle & Fal- 


|garone| ( |2012| ; Padoan et al. (20131 and recent theoretical 
modelling can be found in Federrath, Klessen fc Schmidt] 
(20081; Federrath (20131; Federrath & Klessen (20131. 


Observational evidence to test the simulation results is 


slowly growing. 

Padoan et al. (1997 

1, Goodman, Pineda & 

Schnee (20091, 

Kainulainen et al. 

( 

20091, Schmalzl et al. 

(20101, Schneider et al. (20121 and 

Kainulainen et al.|(|2014|) 


obtained lognormal column density PDFs for the extinction 
through dust in molecular clouds in the Milky Way (MW). 
The latter authors discuss the relationship between the den¬ 


sity PDFs and star formation. Schneider et al. (20131 pre¬ 


sented lognormal PDFs of 12CO(1-0) integrated intensities 
for four dust clouds in the MW, and Hughes et al. (20131 


found lognormal PDFs of the 12CO(1-0) brightness temper- 
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2 E. M. Berkhuijsen and A. Fletcher 


ature and integrated intensities in M51. Also the PDF of 
the H 2 column densities in M33 is lognormal, but with a 
high-density tail (Druard et al.|2014l. Wada, Spaans & Kim 


(20001 showed that the luminosity function of the H i col¬ 
umn density in the Large Magellanic Cloud is lognormal. 


Hill et al. (2007 20081 found lognormal distributions of the 


emission measures, which are proportional to the integral of 
the electron density squared, in the Wisconsin Ha Mapper 
survey (Haffner et al. 20031 perpendicular to the Galactic 
plane at latitudes |fo| > 10°, and along the plane, with dif¬ 
ferent parameters. The emission measures of the extinction- 


corrected Ha emission from the galaxy M33 (Tabatabaei et 


al. 20071 and the volume density of the dust-bearing gas 


near stars within 400 pc from the Sun (Gaustad & Van Bu 
ren||1993 1 have lognormal distributions. Furthermore, Seon 


(2011 20131 obtained lognormal PDFs for the FUV back¬ 
ground emission from the ISRF, and Bowen et al. (20081 
plotted the volume density distribution of the O vi-line emit¬ 
ting gas in the MW, which is approximately lognormal. Re¬ 
cently, Berkhuijsen & Fletcher (2008 20121 showed that in 
the Milky Way the volume density distribution of the dif¬ 
fuse, ionized gas (DIG) at |6| >5°, and of the diffuse atomic 
gas are close to lognormal. They also discussed differences in 
the PDF dispersion observed for the atomic gas at low and 
high latitudes, and of warm and cold gas. Here we present 
the PDFs of ionized and neutral gas column densities in the 
nearby galaxies M31 and M51, which are close to lognormal, 
and discuss their properties. 

Unfortunately, we cannot separate the variation in gas 
density due to turbulence from variations caused by e.g. su¬ 
pernova explosions, HII regions or gravitation. Therefore, 
our PDFs capture all kinds of density structure along the 
lines of sight (LOS), regardless of the origin. On the other 
hand, high-density structures represent only a small fraction 
of the volume of the radio beams at our spatial resolution 
of about 300 pc. We will interpret the properties of the ob¬ 
served PDFs in terms of an effective LOS, Leff, which is the 
fraction of the total LOS occupied by turbulent cells and 
gas clouds. Furthermore, purely for the sake of clarity in the 
text, we will refer only to clouds rather than to clouds or tur¬ 
bulent cells. In reality, we cannot make a firm distinction. 
We shall return to this point at the end of the paper. 

The paper is organized as follows: The data sets used 
are described in Section and we show the derived PDFs 
and discuss their properties in Section]^ Simple models to 
examine the effects of varying resolution, line-of-sight path 
length and the gas filling factor on the form of the PDFs are 
presented in Section]^ We compare the observed properties 
of the PDFs to those expected from simulation studies in 
Section and in Sectionwe summarize our results. 


2 THE DATA 


M31, the spiral galaxy nearest to us, is at a distance of 
780 ± 40kpc (Stanek & Garnavich 19981. The Sb galaxy 
of large angular extent (more than 3° x 1°) and low sur¬ 
face brightness has an inclination of i = 77.5° (Braun|1991 


Ghemin, Garignan fc Foster|2009 1 and so is nearly edge-on. 


The position angle is PA = 37.7° . Our gas density PDFs 
are ba sed on the Westerbork Hi survey of [Brinks fc Sharif 
(19841 at an angular resolution of 24 x 36 arcsec^, corrected 


for missing spacings, the IRAM 12GO(1-0) survey of Nieten 


et al. (20061 at a resolution of 23 arcsec, and the extinction 


corrected Ha data of Devereux et al. (19941, at a resolu¬ 
tion of a few arcsec, produced by [Tabatabaei fc Berkhuijsen] 
(20101. Nieten et al.| ( 20061 converted the GO data to column 
densities of molecular gas using the constant conversion fac¬ 
tor Xco = 1-9 X 10^°molcm“^ (Kkms“^)“^. Before deriv¬ 
ing the PDF we smoothed each map to an angular resolution 
of 45 arcsec, corresponding to 170 pc x 785 pc along major x 
minor axis in the plane of M31. Each PDF was calculated 
from independent data points (at least 1.67x beamwidth 
apart) above 2x rms noise. 

M51 is a bright Sbc galaxy at a distance of 7.6 Mpc 
(Giardullo et al. 20021. It has an angular extent of about 
lOarcmin, a position angle of PA = —10° and an incli¬ 
nation of f = 20° (Tully 1974 Pety et al. 20131. For the 


PDFs we used the Hi data from the THINGS survey (Wal¬ 
ter et al. 20081 at a resolution of 5.82 x 5.56arcsec^, the 


BIMA 12GO(1-0) survey of Heifer et al. (20031 at a resolu¬ 
tion of 5.8 X 5.1 arcsec^, which was corrected for missing 
spacings using single-dish observations, and the Ha map 
of Greenawalt et al. (19981 at a resolution of a few arc¬ 
sec. All maps were reduced to a common field of about 
8.5 X 10.5 arcmin^ covering most of M51 apart from the ex¬ 
tended southern outer arm in Hi. They were then smoothed 
to a resolution of 8 arcsec corresponding to 295 x 315 pc^ 
along major x minor axis in the plane of M51. We converted 
the GO data to column densities of H 2 using the conversion 
factor 4.75 X10^® molcm“^ (K kms“^)“^, one quarter of that 
for the Milky Way ( Guelin et al.|[TOT5 |. Of each map, only 
independent data points above 2x rms noise were used to 
produce the PDF. 


3 PROBABILITY DISTRIBUTION 
FUNCTIONS 

In this Section we present the observed surface density PDFs 
□ for M31 and M51 and characterise them by fitting lognor¬ 
mal distributions to the histograms. Since there are various 
ways to paramaterise a lognormal we shall briefly describe 
the exact equations we use. To allow for easy comparison 
with [Berkhuijsen &: Fletcher| ( |2008P and other work on ob¬ 
served PDFs we use common, or base 10, logarithms. The 
lognormal probability distribution function for the random 
variable x is given by 

A(a;|^, a) = exp [-(logn, x - [if /2a^] dx, (1) 

axy/Zn 

which is equivalent to a Gaussian distribution for logiQ(a:), 
^[loglo(2;)lM,C^l = 

exp [-(login * - M)V2crVlogio x. (2) 

cry zTT 

We use Eq. to derive [i and a from fits to histograms of 
login (a:), taking into account the logarithmic bin-width to 
ensure that the integral of the PDF is unity. It is important 


^ We call our PDFs of the gas densities probability distribution 
functions, rather than probability density functions, to avoid re¬ 
ferring to the “probability density of the density”. 
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Figure 1. Probability distribution functions of the surface den¬ 
sities of ionized and neutral gas in the radial interval 0' < 
R < 60' in M31. Top: Ho (in cm“®pc), middle-. Af(Hl) (in 
atomcm“2), bottom: Af(H 2 ) (in 10^® molcm“®). Dashed 
lines show the lognormal fits given in Table 


Figure 2. Probability distribution functions of the surface 
densities of ionized and neutral gas in M51. Top: Ho (in 
cm“®pc), middle: Af(Hl) (in 10^® atomcm“®), bottom: A1(H2) 
(in 10^® molcm“®). Dashed lines show the lognormal fits given in 
Tabled 


to note that and a are the mean and standard deviation 
of the Gaussian distribution of log]^Q(a:); the mean and stan¬ 
dard deviation of the distribution of x are given by 

Elx] = exp {fi -I- cr^/2), 

and 

S[a;] = \/(expo® — 1) exp (2/i -|- cr®), 

respectively. Thus the variance of x depends on both the 
mean and standard deviation of log]^Q(x). We shall try to 
avoid confusion by calling a the dispersion of the PDF and 
using Xo = 10'' to describe the peak of the PDF (note 
that 10'“ is the median of the distributions of both x and 
logj^g(a;)). For the Gaussian PDFs of log^Q(a:), which is the 
form we show in our figures, small (large) values of a result 
in a narrow (broad) Gaussian PDF. 


3.1 PDFs of surface densities 

In Fig.[2we present the M31 PDFs of the emission measures 
of the ionized gas and of the column densities of H i and H 2 
for the radial range 0' < R < 60' in the plane of the galaxy, 
corresponding to 0 < i? < 13.6 kpc. In Table we list the 
parameters of the lognormal fits together with those for the 


area containing the bright emission ring (35' < R < 60' or 
7.9 < R < 13.6 kpc) and for the weaker emission interior to 
this ’ring’ {O' < R < 35'). For Hq and H 2 the fits of the 
PDFs for the subregions are nearly the same, but for H i the 
position of the peak of the PDF, Xo, for the bright ’ring’ 
is at a significantly higher density than for the inner region 
and the dispersion is smaller. This difference in Xq causes 
the low-density excess — and hence the asymmetry — in 
the total PDF and its larger dispersion compared to that of 
the subregions. 

Figurej^shows the M51 PDFs of the emission measures 
of the ionized gas and of the column densities of H i and H 2 at 
radii R < 3'. The parameters of the lognormal fits are given 
in Table Again the H i PDF has an excess at low densities. 
The peak at high densities is due to gas in the outer spiral 
arms. We have tried to construct separate PDFs for spiral 
arms and the interarm regions by using a spiral-arm mask; 
this was created using the coefficients of a wavelet transform 
of the CO map with a scale of about 1 kpc, the positive 
coefficients tracing a clear two-armed spiral. Although the 
dispersions of these PDFs could not be reliably determined, 
the values of Xo of the arm PDFs are significantly higher 
than of the PDFs for the interarm regions. Recently, |Hughes| 


et al. (20131 found that PDFs of the integrated 12CO(1-0) 
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4 E. M. Berkhuijsen and A. Fletcher 


intensities in M51 for arm regions are clearly wider than 
for interarm regions, and have an about 25 per cent larger 
Xq. They also presented PDFs for the inner 4.5' x 2.8' of 
M51. At a spatial resolution of 53 pc, using all pixels with an 
intensity above 3 times the noise, they obtained a dispersion 
of a — 0.62, the same as our value for a larger area of M51 
observed at a lower resolution. 

The PDFs of Ha, Hi and H 2 for M31 and M51 pre¬ 
sented in Figs. and are close to lognormal. That the 
PDFs of H 2 are lognormal may seem surprising in view of 
the recent work of |Shetty et al. (2011 \. In their simulations 
of molecular clouds, the authors found that the PDFs of 
12CO(l-0)-intensities and of CO column densities are gen¬ 
erally not lognormal and have different shapes when CO is 
optically thick. But at low densities and for the model of a 
Milky Way cloud of medium density the PDFs are lognor¬ 
mal. The box size in their simulations is 20 pc, much smaller 
than the typical resolution of 300 pc used for the PDFs of 
the gas surface densities in M31 and M51. Hence, the line 
intensity seen by the radio beam is the average intensity of a 
large sample of individual clouds. Saturated lines will hardly 
be observed because the filling factor of the densest clouds 
is very small. Therefore the observed CO PDFs become log¬ 
normal and will give a fair representation of the PDF of the 
H 2 column density. 

Comparing Figures and we note the following: 

(a) Both PDFs of Ha show a low-intensity cut-off at the 
sensitivity limit. A similar cut-off occurs in the PDF of H 2 
for M31. 

(b) The Ha PDF for M51 extends to higher emission 
measures than that for M31 and the emission measure of the 
maximum is about 35% higher than that of M31. Since the 
M51 data were not corrected for extinction and the inclina¬ 
tion of M51 is smaller than that of M31, this means that the 
face-on surface brightness in Ha of M51 is much larger than 
that of M31. 

(c) The steep decrease seen in the Hi PDFs at high 
column densities may be attributed to opacity in the Hi 
lines. In scatter plots of star-formation rate against Hi mass 
surface density a cut-off near 10 Mq pc“^ is visible in M31 


(Tabatabaei & Berkhuijsen 

2010|, M51 and many other 

galaxies ([Bigiel et al.|2008|). 

Braun et al.N2009| have shown 


that after correcting the HI column densities of M31 for 
opacity the apparent cut-off vanishes. 

(d) The Hi column densities of M31 extend to higher 
values than those of M51 and the maximum of the PDF 
of M31 occurs at a column density that is 5 times that of 
M51. On the sky, M31 is much brighter in Hi than M51, but 
the face-on value of Xq, defined as Xocosi, of M31 is only 
about 1.2 times that of M51. The Hi PDF for M51 has a 
pronounced low-density tail, while that for M31 only shows 
a mild low-density excess. The low-density tail of the M51 
PDF may be a relic of the last encounter with the companion 
galaxy NGC5195 ( [Howard fc Byrd||1990P that gave rise to 
the extended Hi arm in the south-east ( Rots et al.||1990 K 

(e) The H 2 column densities of M51 are much higher 
than of M31: the face-on value of Xq of M51 is 11 times that 
of M31. The star formation rate of M51 is five times higher 
than that of M31 (Tabatabaei et al. 20131, which may be 
related to the higher column densities reached in M51 than 
in M31 (Kainulainen et al. 20141. H 2 is the dominant gas 
phase in the inner part of M51 (Garcla-Burillo et al.|[T9^ 


M31 35’<R<60’ 



log (VOLUME DENSITY) [at/cm'’] 

Figure 3. Probability distribution functions of mean volume den¬ 
sities along the line of sight (in cm“^) of neutral gas in the radial 
interval 35' < R < 60', the bright emission ’ring’, in M31. Top: 
(uhi), middle: (2nH2), bottom: (n) = (njji) + (2nH2). Dashed 
lines show the lognormal fits given in Table [^ 


Schuster et al. 20071, while in M31 the atomic gas is the 
dominant component (Dame et al.||1993 Nieten et al.|2006|. 


3.2 PDFs of volume densities in M31 

We calculated the mean Hi volume densities, (nni), from 
the column densities using the formula for the Hi scale 
height increasing with radius, h{R), given by [Braun] ( |1991| ), 
scaled to the distance D = 780 kpc. This gives /i(lO') = 
240 pc and h(60') = 425 pc. The volume density then is 
(n-Hi) = A(Hi)/I/ where the line of sight L — 2h{R)/cos i. 
Assuming a scale height of the molecular gas equal to 
h{R)/2 yielded the mean volume densities of the molecu¬ 
lar gas, (n 2 H 2 )- The volume density of all neutral gas then is 
(n) = (jiHi) + (n 2 H 2 )- As an example, we present the PDFs 
of these volume densities for the ’ring’ in Fig. and the 
lognormal fits for all regions in TableThe Hi PDF Is nar¬ 
rower and more symmetric than that for the full area shown 
in Fig. [^because the low-density extension is not included. 
It is shifted by a factor of 2.4 w.r.t. the PDF for the inner 
region. Comparison with Table shows that the dispersions 
of the PDFs of (tihi) and (n 2 H 2 ) are the same as those of 
the column density PDFs within the errors (see Sectionj^for 
a simple explanation for this equivalence.). These estimates 
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Gas Density PDFs for M31 and M51 5 


Table 1. Lognormal fits to the PDFs of gas surface densities in M31, HPBW 45". 


X 

Units 

m 

N 

Position of 

maximum 

Yo 

Dispersion 

(7 


Af(Hi) 

10 ^® atomcm“® 

0-60 

1249 

2.16 ±0.03 

145 ± 10 

0.27 ± 0.02 

12.7 



0-35 

368 

1.86 ±0.01 

72 ±2 

0.23 ±0.01 

0.8 



35-60 

881 

2.32 ±0.02 

209 ± 10 

0.19 ±0.01 

5.5 

Af(H2) 

10^® molcm“® 

0-60 

971 

1.19 ±0.04 

15.5 ± 1.5 

0.41 ± 0.02 

2.9 



0-35 

250 

1.11 ± 0.05 

12.9 ± 1.6 

0.39 ±0.04 

1.5 



35-60 

721 

1.22 ±0.04 

16.6 ± 1.6 

0.40 ±0.03 

3.0 

Ho 

cm“^ pc 

0-60 

1097 

0.88 ±0.03 

7.6 ±0.5 

0.34 ± 0.02 

2.9 



0-35 

397 

0.93 ±0.03 

8.5 ±0.6 

0.34 ± 0.02 

1.2 



35-60 

700 

0.85 ±0.03 

7.1 ± 0.5 

0.34 ± 0.02 

2.3 


N is the total number of data points; the error in each bin Si is estimated as Si = y/Ni for (number of bins—2) degrees of freedom; is 

the reduced chi-squared goodness of fit parameter. 


Table 2. Lognormal fits to the PDFs of gas surface densities in M51 at ij < 3', HPBW 8". 


X 

Units 

N 

Position of 

maximum 

Yo 

Dispersion 

<j 

X^ 

X(Hi) 

10^^ atom cm~^ 

2015 

1.46 ± 0.05 

28.7 ± 3.3 

0.51 ± 0.05 

6.4 

X(H2) 

10 ^® molcm“® 

282 

1.60 ±0.05 

39.4 ± 4.5 

0.64 ±0.05 

1.1 

Ha 

cm“® pc 

728 

1.01 ± 0.04 

10.2 ± 1.0 

0.59 ±0.03 

1.4 


N is the total number of data points; the error in each bin Si is estimated as Si = \/7Vi for (number of bins—2) degrees of freedom; x^ is 

the reduced chi-squared goodness of fit parameter. 


for the mean volume density PDFs are useful in that they 
allow us to compare the properties of M31 with those of the 
Milky Way in Section |5.2| 


4 DEPENDENCE OF DENSITY PDF 
PARAMETERS ON PATHLENGTH, 
RESOLUTION AND FILLING FACTOR 


The maps used to compute the column and average volume 
density histograms presented in Section all have a reso¬ 
lution of a few hundred parsec. The presence of a lognor¬ 
mal density PDF does not in itself indicate that the gas is 
turbulent, the lognormal PDFs are merely compatible with 
this assumption (Tassis et al.|20l'0 Hennebelle & Falgarone 


|2012[ ). One could consider two contributions to the column 
density PDFs: one from the density structure arising from 
compressible turbulence in the ISM and the other from the 
densities of specific, discrete objects that are not caused by 
turbulence, such as SNRs and Hii regions (both of which 
may, however, be drivers of turbulence). How to separate 
these contributions to the PDF is not clear. Both types of 
structure have sizes less than the beam. The dominant en¬ 
ergy input driving turbulence in the ISM is probably due 
to stellar cluster winds and superbubbles produced by clus¬ 
tered supernovae, with a typical scale in the range 50-500 pc 
( Elmegreen fc Scalo|2004 Section 3). Many other drivers of 
turbulence, such as spiral shocks, the magneto-rotational in¬ 
stability and proto-stellar jets have also been proposed (see 


[Federrath fc Klessen|[20T3l for a discussion of these in the 
context of PDFs). The H 2 at the largest scale is collected in 
giant molecular clouds (GMCs), with a typical size of 50 pc 
to 100 pc in M51 


and the Ha has a diffuse component as well as higher den- 


sity H II regions of typical size 30 pc (Habing Israel 

1979 

Scoville et al.|2001 Azimlu, Marcinak & Barmby|2011 

. The 

thickness of the diffuse gas layer is of the order of 500 pc and 


that of the H 2 150 pc, so each telescope-beam cylinder will 
contain one or many turbulent cells or GMCs or Hii re¬ 
gions. Furthermore, the ISM consists of several phases, so 
the emission that we use to construct the average volume 
density PDFs is not uniformly distributed along the total 
LOS but only arises from a fraction of this distance. In order 
to understand how averaging over multiple cells or clouds, 
both along and across the beam, and how the filling factor of 
the ISM occupied by the phase we observe affect the density 
PDFs we constructed a simple model. 


(Blitz et al.|[2007 


Colombo et al.| 2014 


4.1 A simple model for the random gas density 
distribution in the ISM 

Independent random numbers drawn from a lognormal dis¬ 
tribution with parameters p = 0, cr = 1 — see Eqs. 0 
and 0 — were distributed over a Cartesian (a;, y, z) grid of 
(1000 X 1000 X L) mesh-points. The value of each mesh-point 
represents the gas density of a single gas cloud and the unit 
of distance in the model is thus the cloud size. The parame¬ 
ter L represents the line-of-sight and for L > 1 the random 
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Table 3. Lognormal fits to the PDFs of gas volume densities in M31, HPBW 45". 


A 

Units 

R(') 

N 

Position of 

maximum 

-^0 

Dispersion 

cr 


(UHI> 

atom cm~^ 

0-60 

1457 

-0.75 ±0.03 

0.178 ±0.012 

0.25 ± 0.02 

12.5 



0-35 

459 

-1.00 ± 0.02 

0.100 ± 0.005 

0.23 ± 0.01 

1.6 



35-60 

998 

-0.62 ±0.02 

0.240 ± 0.011 

0.19 ± 0.01 

5.7 

("2H2> 

atom cm~^ 

0-60 

920 

-1.05 ±0.03 

0.089 ± 0.006 

0.39 ± 0.02 

2.6 



0-35 

248 

-1.08 ±0.05 

0.083 ±0.010 

0.41 ± 0.03 

1.3 



35-60 

672 

-1.04 ± 0.03 

0.091 ± 0.007 

0.37 ± 0.02 

2.0 

(n> 

atom cm~^ 

0-60 

1451 

-0.61 ± 0.02 

0.245 ±0.012 

0.30 ± 0.01 

5.1 



0-35 

457 

-0.82 ±0.02 

0.151 ± 0.007 

0.32 ± 0.01 

1.1 



35-60 

994 

-0.51 ± 0.01 

0.309 ±0.007 

0.25 ± 0.01 

2.6 


is the reduced chi-squared goodness of fit parameter: the error in each bin 5i is estimated as 5i = y/Nl for (number of bins—2) 

degrees of freedom. 




FWHM 




Figure 4. (a) The effect of an increasing LOS on the dispersion 
(f7 in Eq. 0 ) of the best fit lognormal PDF. The dashed line 
shows the best-fit power-law. (b) The effect of smoothing with a 
2D Gaussian, of unit mean and full-width half-maximum shown 
on the horizontal axis, on a. Circles are for models with L = 1 and 
triangles for L = 10. The inset panel shows the best-fit power- 
laws to the data points where FWHM ^ 1. In both panels the 
modelled column densities were calculated using the simple model 
described in Section [4| 


numbers were summed along the 2 :-axis to give a map of the 
column density of the model, with 1000 x 1000 lines of sight. 
An average volume density along each line of sight was also 
computed by dividing the column density by L. Convolu¬ 
tion of the column density with a 2D Gaussian of a given 


Figure 5. (a) The effect of varying fractional volume, on the 
dispersion {a in Eq. Q) of the best fit lognormal PDF of the 
average volume density. Lines are shown for different LOS. (b) 
(T versus the effective LOS Lgff as defined in the text. The best- 
fit power-law for > 1 is shown. In both panels the modelled 
average volume densities, the column density divided by the LOS, 
were calculated using the simple model described in Section 

full-width half-maximum (FWHM) was used to model the 
consequences of smoothing by a beam of size D = FWHM. 

Making {1 — fv) of the mesh-points empty allowed us to 
control the fractional volume occupied by the clouds, /„ = 
F/(1000x lOOOxL), where V is the number of meshpoints for 
which the gas density is non-zero. The empty mesh-points 
were chosen at random. We shall refer to /„ as the volume 
filling factor but stress that there are at least two other 
quantities that are also called filling factors. Both measure 
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the “dumpiness” of a given gas phase, either averaged over 
the entire volume or only the volume occupied by that phase, 
and calculated using the ratio (n)^/(n^) where n is the gas 
density and (...) represents the appropriate average: note 
that the relation between these filling factors and /„ is not 
straightforward (see Gent et al. 2013 Sect. 5). Under the 
assumption that the gas is confined to clouds which all have 
the same density, which is commonly made when deriving 
filling factors from observational data (see e.g. [Berkhuijsen 


fc Fletcher|2008 1, the fractional volume occupied by the gas 


and the volume filling factor defined as {n)^/{rd) are the 
same. We adopt this assumption in what follows; in other 
words we assume that the fractional volume and the volume 
Hlling factor are one and the same quantity. 

Normalised histograms of the column or average-volume 
density gave the model PDFs, and the best-fitting lognor¬ 
mal function to the PDF, using a least-squares method, was 
used to compute /r and a for each model. Ten independent 
realisations were calculated for each set of model parameters 
and the resulting /r and cr averaged. For L = 1, /„ = 1, and 
perfect resolution (no smoothing), the model PDFs have the 
same parameters, jj, — 0 and u = 1, as the lognormal dis¬ 
tribution used to populate the model (within a numerical 
error due to the finite size of the mesh and the constraints 
of the fitting algorithm). We used the model to explore how 
cr, as this is by far the most useful of the two parameters in 
interpreting observations, varies as L, /„ and D change. 


Figure, with a oc D~^ irrespective of L. This power law be¬ 
haviour means that when D > 4-5 clouds for shallow L and 
D > 2-3 clouds for longer L, the effect on cr of further mod¬ 
erate smoothing is marginal: the resolution of observations 
of nearby galaxies is typically in the range 1 < D < 10. In 
principle, if observations are available with a resolution of 
D < Smax and repeated smoothing results in u oc D~^ be¬ 
yond a certain beamwidth, one could recover Smax from the 
column density PDFs. 


4.5 PDF dispersion and the filling factor 

Figure l^a) shows that as the hlling factor increases, i.e. as 
the gas occupies more of the total volume, the dispersion 
of the htted lognormal PDFs for the average volume den¬ 
sity decreases linearly. However, the effect on cr depends on 
both L and /„. If /„ is a constant for all lines of sight, then 
the power law scaling a oc L°‘ with a = —0.5 still holds. 
However, if /„ depends on L this scaling can change. This 
motivates us to dehne the effective pathlength 

Teff = /ai (3) 

which is the fraction of the total line of sight that is occupied 
by clouds or turbulent cells. Figure shows that a scales as 
a power law a oc for Lets > 1. The best ht power law has 
an index a = —0.5, as in the models where /^ = 1 shown in 
Fig.ga). 


4.2 PDF dispersion and the pathlength 

Figurej^a) shows how the dispersion of the lognormal PDF, 
cr, changes as the pathlength along the line of sight, L, is 
increased. Here /„ = 1 and no smoothing is applied. For L = 
1 the model recovers ct = 1, as required by the parameters 
of the original distribution of densities. As L increases cr 
decreases as a power law cr oc Z/“. The fitted power law in 
Fig. 0 a) has a = —0.4 for the range 1 < L < 64, but as the 
lower end of the range in L increases (i.e. for L ;§> 1) we find 
that a —^ —0.5. This scaling is in agreement with |Fischera| 
& Dopita ( 2004|) who studied the relationship between cr of 


the column density PDF and the ratio of L to the size of a 
turbulent cell, T/smax: for thick discs, i.e. L/sma.^ S> 1, they 
derived cr oc (I//smax)~°'^. In our simple models Smax = 1 
and so L/smax > 1. 


4.3 PDFs of column and average volume densities 


The scaling of a with L for the lognormal PDFs of the aver¬ 
age volume density, dehned as (n) = E/L, follows the same 
power law behaviour as the column density PDFs: for n we 
also find that a oc for L 3> 1. This is not surprising, as 

dividing a lognormal by a constant changes fj, but not cr (e.g. 
Aitchison fc Brown||1957 Theorem 2.1). In these models we 


set /„ = 1 and applied no smoothing. 


4.4 PDF dispersion and beam smoothing 

Figure |^b) shows how cr changes under smoothing of the 
model column (and volume) densities, where we have fixed 
fv = 1- Smoothing reduces cr and for beamwidths D ^ 1 
the scaling follows a power law, shown in the inset to the 


4.6 Summary of model results 

For convenience we summarise the results obtained using 
this simple model for a lognormal distribution of gas den¬ 
sity in the ISM. We shall use these results to interpret the 
observations. 

(i) The dispersions of the lognormal column density and 
average volume density PDFs both scale as a power law with 
the total line of sight, a oc L“. When the hlling factor /« = 1 
both have a = —0.5. 

(ii) Smoothing by the telescope beam, when the beam is 

bigger than the size of a cloud, results in a power law scaling 
cr oc D~^. For D > 5smax in a thin layer (L/smax ~ 1) 
and D > 2smax in a thicker layer ~ 10), further 

smoothing only produces a weak effect on cr. 

(iii) a is linearly related to the hlling factor of the gas 
but the magnitude of the effect also depends on L. Dehning 
the effective pathlength Teff = fvL results in the scaling 
cr « LJ-7 


5 DISCUSSION 

Two diherences between the PDFs shown in Figs. [lj|^and 
Tables m are immediately obvious: the M51 PDFs are 
signihcantly wider than the M31 PDFs and for both galaxies 
the H 2 PDF is wider than the Hi PDF. Which factors could 
cause these diherences? 

Simulations of the ISM have shown that the dispersion 
of a PDF depends on variables like disk mass, temperature, 
magnetic held strength, length of the line of sight and den¬ 
sity scale. Furthermore, we showed in Sectionj^that also the 
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filling factor plays a role. Below we consider each of these 
factors separately. 

The width of the PDF increases with the gas density 
in the disc ( |Tassis|2007| |Wada fc Norman|2007[ ). In the HD 
simulations of Wada & Norman (20071 an increase of a fac¬ 
tor of 10 in disc mass leads to an increase in the dispersion 
of 30%. As the mean face-on column density of the total gas 
in M51 is about 5 times higher than that in M31 (calculated 
from the Hi and H 2 PDFs with (N) = Aq exp(crf„/2) and 
added ( Ostriker et al.|2001 \), this could account for a differ¬ 
ence of about 15 per cent in the dispersions of the gas PDFs 
for the two galaxies. 

In the MHD simulations of ide Avillez fc Breitschwerdtl 
(20051, the dispersion of the PDF of cold gas is wider than 
that of warm or hot gas. This was confirmed by the HD 
simulations of |Robertson &: Kravtsov] ( |2008P and [Gent et al.| 
(20131. Since the molecular gas is generally cooler than the 
atomic gas, the temperature difference could cause the H 2 
PDF to be wider than the Hi PDF. Unfortunately, we can¬ 
not estimate the expected difference in dispersion because 
the temperature ranges in the figures of these simulation 
papers are too large. 


Molina et al. (20121 simulated the development of 


molecular clouds with time using a MHD code. They de¬ 
rived relations between the dispersion of the PDF, the Mach 
number M and plasma /?, the ratio between thermal and 
magnetic pressure. They found that increasing M, i.e. more 
cooler and denser gas, widens the PDFs of the clouds (see 
also I Federrath et ah] |2008[ [Price et al.| |2011[ |Federrath| 


20131, consistent with our observations, whereas decreasing 


/3 makes the PDFs somewhat narrower. We can estimate the 
ratio of mean /3 in M51 to that in M31 by comparing the 
thermal and magnetic pressures in the two galaxies. As the 
total thermal pressure is dominated by the ionized gas, we 
use Xo of the Ha PDFs (Tables 1 and 2) and the lines of 
sight L given in Tables 4 and 5 to calculate the rms elec¬ 
tron density \/(niy for each galaxy. This is 2.9 times higher 
in M51 than in M31 suggesting that in M51 the thermal 
pressure is about three times higher than in M31. As the 
magnetic field in M51 ( Fletcher et al.||2011 | is about three 
times stronger than in M31 ( Fletcher et al.|2004 |, the mag¬ 
netic pressure in M51 is about nine times higher than in 
M31; hence, 13 in M51 is about one third of that in M31. 
Molina et al. (20121 set /3 oc M~^ in their simulations of 


molecular clouds. If this relationship also holds in less dense 
gas phases, a three times smaller (3 would imply a \/3 = 1.7 
times higher mean M in M51 than in M31. Since M affects 


the PDF dispersion more strongly than fi (see Molina et al. 


20121, the PDFs for M51 will be wider than those for M31, 


in agreement with our observations. Separate measurements 
of M and /? will be necessary to see how large this effect is 
compared to other influences on the width of the PDF. 

The scale of the density fluctuations influences the dis¬ 
persion of the PDF because an increase in the number of 
clouds seen by the radio beam leads to a decrease in the dis¬ 
persion as extremes are averaged out (Ostriker et al.pOOl 


Vazquez-Semadeni fc Garcia|200T I. In Section we demon- 


strated that this effect may occur across the beam area as 
well as along the LOS. In Section [5.Ij we show that the an¬ 
gular resolution hardly influences the observed dispersion, 
and in Sectionthat the effective length of the LOS is the 
dominant factor in shaping the observed PDFs. 


M31 0’ < R < 60’ 



Figure 6. Dependence of the dispersion a of the PDFs of the 
column densities Af(Hl) (dots) and A(H 2 ) (triangles) on the half¬ 
power beamwidth (in arcseconds) for the radial interval 0' < iJ < 
60' in M31. A beamwidth of 45" corresponds to 170 pc X 785 pc 
in the plane of the galaxy. The lines are weighted least-squares 
fits to the data points. 


5.1 Does the observed dispersion depend on 
beamwidth? 

In Section we showed that the dispersion a of the col¬ 
umn density PDFs varies strongly under smoothing when 
the smoothing beamwidth is smaller than about 5 clouds, 
but that once the beamwidth is larger than this the effect of 
further smoothing becomes negligible. Thus the behaviour 
of the best fit a under repeated smoothing can, in principle, 
be used to set a rough upper limit on the correlation size of 
the clouds. 

We investigated this point for the column densities of 
Hi and H 2 of M31. To this end we smoothed the Hi and 
H 2 maps with Gaussian functions to angular resolutions of 
24arcsec x 36arcsec (indicated by 30arcsec), 90arcsec and 
180 arcsec, and calculated the PDFs. Including the PDFs at 
45 arcsec and the original resolution of 23 arcsec of the H 2 
map, the range spans a factor of 7.8 in scale. We found that 
for the full region R < 60' as well as for the subregions 
the dispersions of the PDFs are the same within the errors 
for both Hi and H 2 (see Fig. [^. Hence, the dispersion of 
the gas PDFs for M31 is independent of smoothing scale 
between about 90 pc x 400 pc and 700 pc x 3150 pc in the 
plane of M31. 

Since the LOS through the disc of M31 is large (on aver¬ 
age 3000pc for the Hi disc), Fig.|^b) suggests that smooth¬ 
ing will not affect o when the FWHM is wider than about 
2 clouds. We can use this to make an estimate of the upper 
limit to the size of the clouds in M31. The smallest of the 
smoothing scales we consider, 90 pc x 400 pc, must contain 
at least 2^ = 4 clouds, giving an estimate for the upper limit 
of the size of the clouds of Smax < -\/(400 x 90)/4 = 95 pc. 

Note that the linear resolution of the M51 observations 


© RAS, MNRAS 000,[2{^ 



















































Gas Density PDFs for M31 and M51 9 


is about 300pc, similar to that of the best resolved Hi ob¬ 
servations of M31. Therefore, the difference in dispersion of 
the PDFs of the gas surface densities for M31 and M51 can¬ 
not be caused by a difference in angular resolution, provided 
the scales of the clouds are the same. 


5.2 Line-of-sight effects 

In the case of column density PDFs, Figure|^a) shows that 
at a fixed resolution the fitted cr will be smaller by about a 
factor of 2 for a 4-times longer LOS. In particular, we found 
that a decreases with We observe this effect in the 

PDFs of M31 and M51 in several ways. 

First, the PDFs of M31 (Fig. Q have a smaller dis¬ 
persion than those of M51 (Fig. [^. The LOS through the 
strongly inclined galaxy M31 is typically about 4 times 
longer than the LOS through M51 that is nearly seen face- 
on, provided that the scale heights are similar. Assuming 
that the cloud scales in the two galaxies are comparable, we 
then expect more clouds along the LOS through M31 than 
along the LOS through M51, and hence smaller dispersions 
of the PDFs for M31 than for M51. 

Second, for both M31 and M51 the PDFs of H 2 are sig¬ 
nificantly wider than those of Hi. Since the molecular gas is 
more concentrated in the disc mid-plane than H i (due to the 
larger gravitational forces) and molecular clouds are gener¬ 
ally denser and smaller than Hi clouds, the LOS through 
molecular gas are shorter than through Hi gas leading to 
a larger dispersion of the H 2 PDF than of the Hi PDF. 
Also the Mach number of the H 2 gas is higher than of the 
HI gas (because the temperature is much lower whereas the 
turbulent velocity is similar), causing a higher compression 
of the gas in the molecular clouds by turbulent shocks and 
larger dispersions of the cloud PDFs ( [Molina et al.| ( |2012| ) 
and references therein). However, this effect may be reduced 
if the LOS intersects many clouds, or if the beamwidth is 
larger than the typical cloud size: in this case the peak of 
the PDF may move to a higher density but the dispersion 
can be unchanged. Which of these two effects is responsible 
for the observed difference in the PDF dispersions requires 
further study. 

Third, the H i PDF for the inner region of M31 is slightly 
wider (a — 0.23 ± 0.01) than that for the bright ’ring’ 
(cr = 0.19 ±0.01) where the LOS is longer and goes through 
more clouds than in the inner region. This interpretation is 
supported by the higher value of Xq in the ’ring’. The same 
trend is visible in the H 2 PDFs (see Table |^. 

Furthermore, the dispersions of gas PDFs derived in nu¬ 


merical simulations of the ISM (e.g. Vazquez-Semadeni & 


Garcia |2001|[de Avillez fc Breitschwerdt|2005[|Wada fc Nor- 


man|200~ I are usually larger than observed PDFs. Berkhui- 

jsen fc Fletcher| ( |2008| ) attributed the difference to the higher 
densities used in the simulations compared to observed den¬ 
sities. However, the much smaller lines of sight available in 
simulations than in observations may be the main cause of 
the larger dispersions of PDFs derived from simulations (see 
also Fischera fc Dopita||2004 (. 

It is possible to quantify the relationship between the 
dispersion of column density PDFs and LOS. Because the 
scale height of the gas in M31 increases outwards, we can 
test this relationship using the PDFs of the average volume 
densities of the sub-regions at different radii given in Table]^ 



Figure 7. Dependence of the dispersion cr of the PDFs of mean 
gas volume densities in M31, the Milky Way and M51 on line of 
sight L (in pc). We assumed an uncertainty of 25% in L. The 
dashed line is the weighted least-squares fit given by Eq. (The 
figure is shown in colour in the electronic version of the paper.) 


In Fig. we show the dispersion of the PDFs of average 
volume density of H i and 2 H 2 as function of the mean LOS L 
of the radial ranges Q' < R < 35' and 35' < R < 60', where 
L = 2h{R )/cost and h{R) is the scale height as derived in 
Section 13.21 

We have added 3 points for the mean volume density 
of diffuse Hi (cool and warm Hi) towards stars in the so¬ 
lar neighbourhood (SN) and for the mean volume density 
of free electrons in the diffuse ionized gas (DIG) at Galactic 
latitudes |6| >5° within a few kpc from the sun as well as 
two points of mean electron density squared obtained from 
emission measures towards pulsars, all taken from jBerkhui-j 
jsen fc Fletcherj (|2008| 20121. We also show the dispersion 


of the Ovi PDF calculated along the lines of sight to 190 
stars using fig. 9 (left panels) in Bowen et al. (20081. 

For M51 scale heights of molecular and ionized gas are 
available. Pety et ^ (2013 \ combined CO surveys observed 
with Plateau de Bure and the 30-m IRAM telescope. An ex¬ 
tensive analysis showed that the molecular gas not only has a 
thin disc with a Gaussian scale height of h — 40 pc, but also 
a thick disc with a scale height of about 250 pc. [^Taking the 
about 8 times higher volume density of the thin disc in the 
mid-plane into account, their sum can be desribed by a single 
exponential with a scale height of about 75 pc. Seonj ([2009 ( 
modelled the Ha emission from M51 observed by |Thilker,| 


Braun fc Walterbos (20001 and derived a thin disc of H ii- 


regions with exponential scale height of 300 pc and a thick 
disc of DIG with 1 kpc scale height. In the mid-plane the 
emission from the thin disc is about 7 times stronger than 
that from the DIG. The sum of the two distributions can 
be approximated by a single exponential disc with a scale 
height of 350 pc. We have added the dispersions of A^(H 2 ) 
and Ha given in Table using the single scale heights ob- 


^ We use the upper value of the range 190-250 pc given by jPetyj 
jet al.j | |2013| l because it applies to the larger radii covered by our 
data. 
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tained above for the LOS L. Details on the data points in 
Fig. are given in Table 

Apart from the A^(H 2 )-point of M51 (that we discuss 
below), the dispersions in Fig. [^systematically decrease with 
increasing L. A weighted fit through the 12 available points 
yields the power law 


log[o-(L)] = (2.0 ± 1.7) - (0.76 ± 0.06) log[L] (4) 
with = 0.3 and L in pc. 

We conclude that the length of the line of sight through 
the gas in M31, the MW and M51 is the main influence 
on the width of the PDFs of average volume density. This 
implies that the number of clouds along the LOS determines 
the dispersion of the PDF. 

Since the data points for the Milky Way and Ha for M51 
in Fig.[^are consistent with the M31 points, the variation in 
cloud sizes must be similar in these galaxies. This conclusion 
is supported by the work of fAzimlu, Marcinak fc Barmby] 
(20111, who found that the size distributions of Hii regions 
in different galaxies are very similar, but the total number 
of H II regions differs depending on the star formation rate. 

We note that most of the data points in Fig. [^forming 
the a - L relation refer to diffuse gas. Much of the Hi in 
M31 is known to be diffuse ( Braun fc Walterbos||199^ and 
the small amount of molecular gas in M31 may ’look’ diffuse 
in our large radio beam because it does not contain many 
small, dense H 2 clouds. 


5.3 Importance of volume filling factor 


The power-law dependence of cr on L presented in Fig. □ h as 
an exponent of —0.76 ± 0.06, which is —0.26 smaller than 
expected from the models In Sect. [^ However, the density 
structure in the ISM contains holes and voids in and be¬ 
tween gas clouds (e.g. Elmegre^|1997| 1998| de Avillez fc 


Breitschwerdt|2005 1 and so we should also take into account 


the volume filling factor, /„, of the gas using Eg. ([3|) . 

Estimates of /„ from observations (Dickey 1993 


Berkhuijsen|1999 1 and simulations (de Avillez fc Breitschw- 


erdt||2005| |Gent et al.||2013[) have shown that Hi gas has 


a larger filling factor than the denser molecular gas. So in 
Fig-lZl/^ will increase with increasing L and LbH will increase 
faster than L. 

For the (ne)-point for the MW at L = 2480 pc in Fig.j^ 


Berkhuijsen fc Muller (20081 derived a volume filling factor 


of fv = 0.08 ± 0.02. Hence the effective LOS producing the 
dispersion of 0.22 ± 0.02 is LeS = (0.08 ± 0.02) * 2480 « 
200 ± 50 pc. Assuming that the relation a oc L~g '^ found 
in Sect. [^ holds in the ISM, we can estimate Leff for the 
other points in Fig. [^ using the observed a in (t/( 0.22 ± 
0.02) = (Leff/(200±50))“° ®. The volume filling factor then 

is fv ~ Leff /L. 


We list the observed and derived properties of the 13 
points in Fig. in Table [^ The errors given are random er¬ 
rors only. The error of 25 per cent in Lgir of the calibration 
point is systematic and affects all LeS and fv in the same 
way. As expected, the values of LeS increase with increas¬ 
ing L varying between about 25 pc in M51 and more than 
200 pc in the nearly edge-on galaxy M31. The filling factors, 
though, vary much less: in M31 between 0.04 ± 0.01 for H 2 
and 0.08±0.02 for the Hi point at L = 3560pc. Apart from 
the deviating H 2 point for M51, also the other filling factors 


are in this range. These values seem amazingly small, but 
they compare well with other information. For the cold (H 2 


and HI) gas at temperatures T < 200K, Breitschwerdt fc de 
Avillez| ( 20051 found from their MHD simulations fv = 0.06. 


Gent et al. (20131 derived from their simulations fv = 0.03 


for cold gas at temperatures T < 500K. In their OVI-line 
study, Bowen et al. (20081 estimate that the size of the gas 
regions causing these lines is about 100-200'pc. Since the 
mean distance is L = 2040 pc, /„ = 0.049 — 0.098 for these 
regions, consistent with our value of /„ = 0.05 ± 0.02. We 
note that at L = 2480 pc a of < rie > is larger than a 
of < rie >, as one would expect. The larger it of < > 

leads to values of Leff and fv that are nearly half of those 
of < Tie >, in agreement with the smaller scale height and 
larger dumpiness of HII regions compared to the DIG. 

In Fig. the H 2 point for M51 is located at a much 
smaller LOS than is expected from the relation formed by 
the other 12 data points. This may be due to the large filling 
factor of fv = 0.15T0.05, which is three to four times larger 
than fv of the H 2 points for M31 (see Table [^. In the inner 
parts of M51 {R < 3') the dominant gas component is H 2 


(e.g. Pety et al. 

20131), while in M31 Hi is the dominant 

component (e.g. 

Nieten et al.||2006l. At the end of Sect. 3.1 


we showed that the mean face-on column density of H 2 in 
M51 is 11 times higher than that of M31, which explains the 
higher value of fv in M51. When the filling factor is high, 
already a small effective LOS contains so many clouds that 
the PDF becomes lognormal. As Lea in M51 is only about 
25 pc, the maximum size of the clouds along the LOS must 
be much less than 25 pc. 

Since the dispersion of the column density PDF is the 
same as that of the PDF of the average volume density (see 
Sect. 1^, we can estimate L of Ha in M31 and of Hi in 
M51 from Eq. [fusing the observed dispersions. This yields 
L = 1750 ± 1500 pc for Ha in M31, similar to that of H 2 
at R > 35 '. Using the same scaling for Lea as above, we 
find Lea = 85 ± 20 pc and fv ~ 0.05 ± 0.02 (see Table [^. 
The value of L implies a scale height of 190 ± 160 pc, the 
same as for the H 2 in the emission ring. This scale height of 
Ha is consistent with the observation of lWalterbos fc BraunI 
(19941 that the z-extent of the DIG in M31 is less than about 
500pc. If Eq.[^and the scaling for Lea also hold for the Hi 
gas In M51, then the observed dispersion in Tablej^indicates 
a LOS of about 1 kpc (see Table |^. The corresponding Lea 
and filling factor of Leii = 37 ± 10 pc and fv = 0.04 ± 0.01, 
respectively, are smaller than in M31, which may be due to 
the small Hi content in M51. 

Interestingly, Seon (20091 presented separate Ha PDFs 
for the DIG and the Hii regions in M51 based on the data 
of Thilker, Braun & Walterbos (20001. Although the scale 
heights differ by a factor of three (300 pc and 1000 pc for 
Hii regions and the DIG, repsectively), both PDFs have a 
dispersion of 0.36 ± 0.03, but the position of their maxima 
differs by about a factor of 7. The LOS through the DIG 
would be 2130 pc and the dispersion of 0.36 ± 0.03 is not 
far from the value of 0.29 ± 0.09 expected from Eq. [^ but 
the point for the Hii regions PDF with an LOS of 640 pc 
disagrees with the value of a = 0.73 ± 0.14 expected from 
Eq. (0. The same dispersion for the two layers indicates 
that Leii is the same and the difference in LOS then means 
that the volume filling factors and/or the maximum cloud 
scales differ. Using the same scaling as before, we find that 
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Table 4. Effective lines of sight Lgff and filling factors 


Gal. 

Gas 

Zone 

N 

L 

Teff |pc] 

fv 

Dispersion 

Ref. 

M31 

2 H 2 

R < 35' 

248 

1370 

60 ± 15 

0.04 ± 0.02 

0.41 ± 0.03 

1 



R > 35' 

672 

1770 

70 ± 15 

0.04 ± 0.01 

0.37 ± 0.02 

1 


Hi 

R < 35' 

459 

2730 

180 ± 40 

0.07 ± 0.02 

0.23 ±0.01 

1 



R > 35' 

998 

3560 

270 ± 60 

0.08 ±0.02 

0.19 ±0.01 

1 

MW 

Hic 

|6| > 5° 

64 

2050 

130 ± 40 

0.06 ±0.03 

0.27 ± 0.03 

3 


H 

\b\ > 5° 

98 

2300 

90± 25 

0.04 ± 0.01 

0.33 ± 0.03 

2 



|6| < 5° 

42 

4630 

300 ± 65 

0.06 ±0.02 

0.18 ±0.01 

2 


A^e 

|6| > 5° 

34 

2480 

200 ± 50* 

0.08 ±0.02 

0.22 ± 0.02* 

2 


Ovi 

| 2 | > 200pc 

: 190 

2040 

100 ± 25 

0.05 ±0.02 

0.31 ± 0.02 

4 



\b\ > 5° 

34 

2480 

115 ± 25 

0.05 ±0.02 

0.29 ±0.02 

2 




157 

1700 

85± 20 

0.05 ±0.02 

0.34 ± 0.03 

2 

M51 

H 2 

R<3' 

282 

160 

24 ± 6 

0.15 ±0.05 

0.64 ± 0.05 

1 


Hq 


728 

750 

28 ±6 

0.04 ± 0.01 

0.59 ±0.03 

5 


N=number of LOS; assumed errors in L: 25%, * a = 0.22 it 0.02 and Left = 200 it 50pc calibrate Left (see text). References: 1. This 
work; 2. Berkhuijsen & Fletcher (2008); 3. Berkhuijsen & Fletcher (2012); 4. Obtained from data in Bowen et al. (2008); 5. Seon (2009). 


Table 5. Estimates of LOS from dispersion of column density PDF. 


Gal. 

Comp. 

L 

[pc] 

Left 

[pc] 

fv 

Dispersion 

M31 

Ha 

1750 ± 1500 

85 ±20 

0.05 ± 0.04 

0.34 ± 0.02 


R < 60' 





M51 

Hi 

1030 ± 900 

37 ± 10 

0.04 ± 0.03 

0.51 ± 0.05 


cr = 0.36 ± 0.03 corresponds to LeS = 75 ± 20 pc giving 
/„ = 0.04 ± 0.01 for the DIG and /„ = 0.12 ± 0.04 for the 
disc of Hii regions. The filling factor of the Hii regions is 
similar to that of the H 2 clouds in M51 suggesting that their 
density structure is similar. 

It is worth noting that as the PDFs for the disc and 
the DIG in M51 are strongly shifted in log(Xo), their com¬ 
bination yields a much wider PDF than that of the separate 
components. Thus our Ha PDF for M51 is not representa¬ 
tive for small regions or specific components. This poses a 
general problem for the interpretation of PDFs for large re¬ 
gions in galaxies containing diffuse, low-density regions as 
well as dense spiral arms. 


5.4 Estimating Mach numbers using PDFs: a note 
of caution 

A connection between the dispersion of the density PDF 
(7 and the Mach number M of interstellar turbulence has 


been identified theoretically (e.g. Padoan et al. 

19971 Passot 

& Vazquez-Semadeni|1998||Ostriker et al.|2001| 

Federrath et 

0 

0 

00 

2010 

Price et al.|2011 

Molina et al.|2012 1 and used 


to make an estimate of M from observations (e.g. Padoan 
et al. 1997 Hill et ar||2008 Berkhuijsen fc Fletcher||2008 1 . 
However, since the magnitude of a depends on the length of 
the LOS L, the telescope resolution D and the volume filling 
factor of the gas /„, caution should be exercised in apply¬ 


ing the theoretical relation to real data. To make a reliable 
estimate of M one needs to know something about the dis¬ 
tribution of the observed gas and so correct the observed a 
for the effects oi L, D and /„. 


5.5 M31: PDFs of dust emission and extinction 


In view of the generally lognormal shapes of the gas PDFs, 
one may ask whether other constituents of the ISM also have 
lognormal PDFs. The dust emission is a good candidate since 
dust optical depth and gas column density are well corre¬ 


lated in the solar neighbourhood (Bohlin, Savage & Drake 
|1978P and in M31 ( |Tabatabaei fc Berkhuijsen||2010p , and 
lognormal PDFs of molecular gas column density and of ex¬ 
tinction through dust clouds in the MW have been reported 


(e.g. Padoan et al. 

1997 Goodman, Pineda & Schnee|2009 

Kainulainen et al.|| 

2009 Schneider et al. 12013 1. 


We calculated the PDFs of the dust emission from M31 
at 24 ^m, 70 /rm and 160 /rm using the Spitzer MIPS maps 
of iGordon et al. (20061 smoothed to a resolution of 45 " 


by Tabatabaei & Berkhuijsen (20101. The lognormal fits are 
given in Table At 24 ^m and 70 /rm the fits for the 3 re¬ 
gions closely agree, as is also the case for the PDFs of N(H 2 ) 
and Ha (see Table [^. The dispersion of 0.34 ± 0.01 of the 
24 fim PDFs is the same as that of Ha and the dispersion 
of 70/rm PDFs agrees with that of A^(H 2 ). As both the ion¬ 
ized gas and the warm dust are mainly heated by young OB 
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Table 6. Lognormal fits to the PDFs of dust emission from M31, HPBW 45”. The fitted function is Y = (V^cr) ^ exp[—(logj^Q X — 
/r)2/2a2]. 


X 

Units 

H(') 

N 

Position of 

maximum 

Ao 

Dispersion 

cr 


24/im 

MJy/sr 

0-60 

1508 

-0.47 ± 0.01 

0.34 ±0.01 

0.34 ±0.01 

1.2 



0-35 

525 

-0.48 ±0.02 

0.33 ±0.01 

0.33 ±0.01 

1.4 



35-60 

983 

-0.46 ±0.01 

0.35 ± 0.01 

0.35 ±0.01 

0.8 

70/rm 

MJy/sr 

0-60 

1426 

0.43 ±0.02 

2.69 ±0.13 

0.39 ±0.01 

1.6 



0-35 

502 

0.41 ± 0.03 

2.57 ± 0.18 

0.37 ± 0.02 

1.8 



35-60 

924 

0.43 ±0.02 

2.69 ±0.13 

0.39 ±0.01 

1.3 

160/rm 

MJy/sr 

0-60 

1516 

1.36 ±0.01 

23.4 ± 0.6 

0.26 ±0.01 

1.8 



0-35 

535 

1.30 ±0.01 

20.0 ±0.5 

0.26 ±0.01 

1.5 



35-60 

981 

1.40 ±0.01 

25.1 ± 0.6 

0.26 ±0.01 

1.4 

Av 

mag. 

0-60 

1406 

-0.04 ± 0.02 

0.91 ± 0.04 

0.20 ±0.01 

6.6 



0-35 

495 

-0.16 ±0.02 

0.69 ±0.03 

0.25 ±0.01 

1.9 



35-60 

910 

0.01 ± 0.01 

1.02 ± 0.03 

0.18 ±0.01 

4.4 


Y is the fraction of sightlines in each bin divided by the logarithmic bin- width d(logjQ X). is the reduced chi-squared goodness of fit 
parameter, with the error in each bin Si estimated as Si = \/Ni and for (number of bins—2) degrees of freedom. 


stars, the agreement in dispersion suggests that the density 
structure of the ionized gas is similar to that of the warm 
dust particles. This is consistent with the good pixel-to-pixel 
correlation between these dust emissions and that of ionized 


gas reported by Tabatabaei & Berkhuijsen (2010 b 


Although at 160 fim the dispersion of the PDFs is also 
the same for each region, the value of Xq for the ’ring’ is 
about 25 per cent higher than for the inner region. The dis¬ 
persion of 0.27 ± 0.01 is close to that of the Hi PDFs (see 
Table [^. Since the dust emitting at 160 ^m is mainly heated 
by the smooth ISRF (Xu & HelouT996l, this suggests that 


the density structures of the cool dust and the Hi gas are 


similar. Indeed, Tabatabaei & Berkhuijsen (20101 found that 


the wavelet correlation on scales below about 1 kpc between 
160 btm and H i is much better than between the other wave¬ 
lengths and Hi. 

We also calculated the PDFs of the dust extinction 
Av through M31 using the map of optical depth at the 
wavelength Ha, thq, derived by |Tabatabaei fc Berkhuijsen| 
( |2010[ ), and Av = 1.234rHa. The parameters of the lognor¬ 
mal fits are given at the end of Table The dispersions 
in the Av PDFs for the subregions closely agree with those 
of the column densities of Hi (see Table [^, the dominant 
gas phase in M31. This indicates that in spite of the radial 
increase in the atomic gas- to-dust ratio the density struc¬ 
ture in the atomic gas and the cool dust are very similar, 
in agreement with the nearly linear correlation between ex¬ 


tinction and Af(Hi) obtained by Tabatabaei & Berkhuijsen 

( [M0| ). 


5.6 Comparison with the Milky Way 

The Hq PDFs for M31 and M51 may be compared to those 


obtained for the MW by Hill et al. (20081 from the Wiscon¬ 
sin Ha-Mapper survey ( Haffner et al.|20'^ b For the PDF of 
the emission measure perpendicular to the Galactic plane, 
EM sin|6|, of the diffuse ionized gas (DIG) at |6| > 10° 


they found a dispersion of 0.190 ± 0.001, which is nearly 
half the value for the DIG in M51 ( Seon|20'09 1 and for M31 
(see Table [l|. H owever, in order to derive the PDF of DIG 
only. Hill et al. (20081 removed all sight-lines towards clas¬ 


sical Hii regions from their sample, whereas in our PDFs 
for M31 and M51 these regions are included. The MW PDF 
including classical H ii regions has a much larger dispersion 
than 0.19 and shows a high-density excess like the Ha PDFs 


for M31 and M51. High-density tails are expected ( Ostriker 
et aI.|200T I if the temperature in the medium decreases with 


increasing density, as is the case in the ionized medium in 
the MW ( Madsen et al.pOOO b 


Hill et al. (2008 \ showed that the dispersion of the PDF 
of EM sin |6| decreases with increasing latitude. |Seon ](S 
found the same effect for the diffuse background emission 
from the MW in FUV, which mainly consists of starlight 
scattered by interstellar dust. Because the LOS decreases 
with increasing |6|, we would expect the opposite behaviour 


(see Sect. 5.21. However, the decrease in a is consistent 
with the observation of [Berkhuijsen fc Muller] ( |2008| ) that 
the spread in (rie) decreases with increasing distance to the 
plane \z\ (also visible in Fig. 5 of Gaensler et al. (20081), 
which is accompanied by an increase in the size of ionized 


clouds. Savage & Wakker (20091 arrived at a similar conclu¬ 


sion based on O vi-line observations. We show below that 
the decrease in dispersion with increasing |&| is due to an 
increase in the volume filling factors of the ionized gas and 
the dust grains, leading to an increase in the effective LOS. 

Hill et al. ( [2008 1 and Seon ( 2011| derived the dispersion 
for the latitude intervals 10° < |6| < 30°, 30° < |6| < 60° 
and 60° < |6| < 90°. From the dispersion and the scal¬ 


ing explained in Sect. |5.3[ we calculated Leir for each case, 
which increases when a decreases. Assuming a maximum 
LOS equal to the scale height of the DIG of L = 1 kpc 


(Berkhuijsen & Muller 2008 Schnitzeler 20121, we obtain 
the volume filling factor /„ = Les/L. Table [(^ shows that 
fv of both the DIG and the scattering dust increases with 
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Table 7. Effective LOS and volume filling factors towards \z\ = 
Ikpc in the MW. 


Data 

\b\ 

(degree) 

cr 

Lefi 

[pc] 

fv 

EM sin|b| 

10-30 

0.19 

260 ± 50 

0.26 ±0.05 


30-60 

0.16 

390 ± 70 

0.39 ±0.07 


60-90 

0.14 

490 ± 90 

0.49 ± 0.09 

Ipuv sin|b| 

10-30 

0.28 

120 ± 20 

0.12 ±0.02 


30-60 

0.23 

180 ± 30 

0.18 ±0.03 


60-90 

0.14 

490 ± 90 

0.49 ± 0.09 


Only random errors due to the error in cr of (ris) are given (see|^. 
Errors in a of EM sin |6| are negligible. Errors in a of Ipjjv 1^1 
are not available. 


increasing |6|. The largest part of the DIG and dust layers is 
sampled by LOS in the interval 10° < |6| < 30°, whereas at 
the highest latitudes mainly DIG and dust at high \z\ near 
the Sun are observed. Therefore, the value of /„ = 0.26±0.05 
for the DIG at 10° < |fo| < 30° best represents the mean 
hlling factor of the DIG seen towards L — 1 kpc perpendic¬ 
ular to the Galactic plane. This value agrees well with the 
value of fv = 0.22 ± 0.05 obtained by [Berkhuijsen, Mitra fc] 


Muller (20061 from emission measures and dispersion mea¬ 


sures towards 157 pulsars at |6| > 10°. These autors found 
an increase of the local Hlling factor fvG) with increasing \z\, 
which was confirmed by Gaensler et al. ( |2008 1 using DM and 
EM towards 51 pulsars with known distance. An increase of 
fv{z) with 1 21 is also indicated by the increase of the mean 
/„ towards higher latitudes in Table 

The dust filling factors for the two lower intervals of |ti| 
are about half of those of the DIG, but at high |b| the values 
of fv are the same. Near the plane the dumpiness of the 
dust seems larger than that of the DIG. On the other hand, 
the values of /„ will be too low if the scale height of the 
diffuse dust is smaller than 1 kpc. 


6 SUMMARY 


We have presented probability distribution functions (PDFs) 
of the surface densities of neutral and ionized gas at radii 
J? < SMn M51 and in the radial ranges 0' < R < 35^ 
35' < R < 60' and 0' < i? < 60' in M31 (Section |3.1 


Figs. [^&[^. For M31 we also presented the PDFs of the 
average volume densities along the LOS of H i, H 2 and total 
gas (Fig. § and derived the PDFs of the dust emission and 
extinction (Table [^. The main conclusions from Sect. 
Sect. and the discussion in Sect. may be summarized as 
follows. 

f. All PDFs are close to lognormal, but their dispersions 
differ (Tables [lj|^. We have investigated which factors de¬ 
termine the dispersion. 

2. The dispersion of the lognormal PDFs of column and 
average volume densities depend on the pathlength L, with 
(7 oc and telescope resolution D, with a oc D~^. The 

PDFs of column and average volume density are also affected 
by the filling factor of the gas; if an effective LOS is defined 
as Leg = fvL then cr oc 


3. The length of the line of sight through the medium, L, 
is the dominant factor shaping the PDF. Short LOS through 
the nearly face-on galaxy M51 cause wider PDFs than long 
LOS through the nearly edge-on galaxy M31, and for both 
galaxies the H 2 PDFs are wider than the H i PDFs because 
the scale height of molecular gas is only half that of atomic 
gas. The contribution of the higher Mach number in the H 2 
to the widening of the observed PDFs needs further study. 

4. The dispersions of the PDFs of Hi and H 2 column 
densities in M31 are independent of the beamwidth for angu¬ 
lar resolutions of 23arcsec to ISOarcsec (Fig.j^ correspond¬ 
ing to linear resolutions in the plane of M31 of 90 pc x 400 pc 
to 700 pc X 3150 pc along major x minor axis. This suggests 
that the maximum cloud size Smax < 95 pc. 

5. The dispersions of the volume-density PDFs of Hi 
and H 2 for M31 (Fig.[^, of diffuse H i and diffuse ionized gas 
for the solar neighbourhood (Berkhuijsen fc Fletcher||2008 


20121, and of diffuse ionized gas for M51 vary with L as ct oc 
6 ±o o6 (Pig indicating that the density structures 
are very similar in M31, M51 and the Milky Way. 

6 . The exponent of the power law a-L relation is steeper 


than the —0.5 expected from analytic calculations (Fischera 
& Dopita 20041 and the models in Sect. We show that 


the difference can be explained by an increase in the volume 
filling factor between dense and less dense regions causing 
an increase in the effective LOS, Leii = fvL. An increase in 
fv also explains the decrease in a with increasing latitude of 
the PDFs observed for the DIG and the scattering dust in 
the Milky Way ( |Hill et al.|2008 Seon|20fT |. 


7. The dispersions of the Av PDFs for M31 closely agree 
with those of the column densities of Hi, the dominant gas 
phase in M31, suggesting that the density structures in the 
cool dust and in Hi are similar. 

8 . The lognormal PDFs we obtained for M31 and M51 
do not, taken in isolation, show that compressible turbulence 
is shaping the gas density distribution. However, they are 
compatible with this interpretation: both in terms of their 
shape and in how their properties vary with beamwidth and 
LOS. By combining the PDF parameters with independent 
measures of, for example, the Mach number, filling factor 
and plasma beta, a deeper understanding of turbulence in a 
multi-phase interstellar medium should become possible. 
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